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Nonlinear gravitational effects of large-scale local inhomogeneities on Cosmic Microwave Back- 
ground (CMB) anisotropies are studied, based on the relativistic second-order theory of pertur- 
bations in nonzero-A flat cosmological models, which has been analytically derived by the present 
author, and on the second-order formula of CMB anisotropies derived by Mollerach and Matar- 
rese. In this paper we derive the components of the CMB anisotropy power spectra in the range of 
I = 1 — 22 which are caused by asymmetric local inhomogeneities on scales of 300 Mpc. Using our 
results it is found that there is a possibility to explain the small north-south asymmetry of CMB 
anisotropies which has recently been observed. 
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I. INTRODUCTION 



In most studies of Cosmic Microwave Background (CMB) anisotropies, the comparison between observed and 
theoretical quantities have so far been done, assuming the linear approximation for cosmological perturbations. 
It seems to be successful enough to determine the cosmological parameters [1, 2, 3]. The present state of 
our universe is, however, locally complicated and associated with nonlinear behavior on various scales, and 
' so the observed quantities of CMB anisotropies may include some small effects caused by large-scale local 
I/"") , inhomogeneities through nonlinear process. Recently it has been reported by Eriksen et al.[4, 5], Hansen et 
t-H ■ al.[6, 7, 8], Vielva et al.[9] and Park[10] that there is a non-trivial north-south asymmetry in various quantities 
*n | about CMB anisotropies. These observational results may suggest the existence of the above effect. 

In this paper we study these nonlinear effects of large-scale local inhomogeneities on CMB anisotropies, based 
on the relativistic second-order theory of cosmological perturbations, which we have recently derived [11] and on 
Mollerach and Matarrese's second-order formula of CMB anisotropies [12]. In Sec. II, we show the second-order 
perturbations in nonzero-A flat cosmological models and the corresponding CMB anisotropies. In Sec. Ill, 
^p" 1 , we derive the expressions for the second-order power (AC;) of CMB anisotropies, assuming a dipole form of 
local inhomogeneities. In Sec. IV, we derive numerically the first-order and second-order anisotropy power 
spectra, and consider the condition that the local inhomogeneities can cause the observed asymmetry of the 
CMB anisotropy, assuming four simple model types of the radial dependence in local inhomogeneities. It is 



found that there is a possibility that the observed decrease of low multipoles of CMB anisotropy[13, 14] also 
may be explained together with the above asymmetry. The derivation of main equations in Sec. II and III are 
shown in Appendix. Concluding remarks follow in Section V. 

II. SECOND-ORDER PERTURBATIONS AND CMB ANISOTROPIES 



First we review the background spacetime and the perturbations which were derived in the previous paper. 
The background flat model with dust matter is expressed as 

ds 2 = g^dx^dx" = a 2 (r])[-drf + <%da;W], (2.1) 

where the Greek and Latin letters denote 0, 1, 2, 3 and 1, 2, 3, respectively, and Sij (= <5* = S^) are the Kronecker 
delta. The conformal time i](= x°) is related to the cosmic time t by dt = a(r])dr]. The matter density p and 
the scale factor a have the relations 

pa 2 = 3(a'/a) 2 - Aa 2 , and pa 3 = p , (2.2) 

where a prime denotes d/drj, A is the cosmological constant, po is an integration constant and the units 
8ttG = c = 1 are used. 

The first-order and second-order metric perturbations Sig^v(= h^) and 52 9p,v{= ^^v)-, respectively, were 
derived explicitly by imposing the synchronous coordinate condition: 



^oo = hoi = and £ 00 = 4i = 0. 



(2.3) 
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Here we show their expressions only in the growing mode: 



h j . = P(j])Fi , 

l\ = P{ V )Li+I#(ri)M>+Q( V )Nl>+Ci, 



(2.4) 



where F is an arbitrary potential function of spatial coordinates x 1 , x 2 and x 3 , h\= S^hu, N, 
and P(rj) and Q(jj) satisfy 

2a' 

P " + f^p'_i = , 

a 

Q" + —Q' = P- b -{P'?- 

CI A 



(2.5) 



The three-dimensional covariant derivative \i are defined in the space with metric dl 2 = S i jdx l dx : > and their 
suffices are raised and lowered by use of Si j . The functions L\ and M- are defined by 



.) ~3F,iF,j — 2.F • F^j + -SijFjFj 



= U 

28 I 



UFuAF - 3<5, 



Fk,F 



kl-f,kl 



(AF) 2 ] } 



and N is defined by 



(2.6) 
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The last term C\ satisfies the wave equation 
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(AF) 2 



UCi = ^{P/a?Gi+ l -[P-l{P>? 



(2.7) 



(2.8) 



where G\ and G\ are second-order traceless and transverse functions of spatial coordinates, and the operator 
□ is defined by 



□0 = g^cp.^ = -a- 2 (d 2 /d V 2 + ^d/d V - A) 



(2.9) 



So C\ represents the second-order gravitational waves caused by the first-order density perturbations. 

The velocity perturbations Si w M and S2 « M vanish, i.e. Si u° = Si u l = and S2 u° = S2 u l = 0, and the 
density perturbations are 



SpIp = -^(-p'-i)af, 

1 pa z V a / 



sp/p 



1 r f 



{^(1 - ^P')(3F,Fj+8FAF) + ^P[(AF) 2 + F M F % 



+ \ [(P 1 ? [(AF) 2 - F M F M ] - ~PP'[4F M F M + 3(AF) 2 ]}. 



(2.10) 



Next let us consider the CMB temperature T = T (0) (l + ST/T), which T (0) is the background temperature 
and ST/T{= SiT/T + S2T/T) is the perturbations. The present temperature is related to the emitted 
background temperature Tj ^ at the decoupling epoch by tJ ^ = (f + z e )To°\ the temperature perturbation 
r = (ST/T) e at the decoupling epoch is determined by the physical state before that epoch , and the present 
temperature perturbations (5T/T) Q is related to (ST/T) e by the gravitational perturbations along the light 
ray from the epoch to the present epoch. The light ray is described using the background wave vector k^ (= 
dx^/dX), where A is the afiine parameter, and its component is k 1 - ^ = (1,— e*), and the ray is given by 
x (°)Ai — (;\ Q _ A)e 4 ], where e l is the directional unit vector. 

The first-order temperature perturbation is 



ST/T = t+- 



dAP'(»7)FyeV. 



(2.11) 
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Using the relation dP/dX = P' and dF/dX — —F yi e l , this equation is expressed as 



where 



ST/T = 6i + e 2 (2.12) 



©i = T-^[(P'F, i ) e -(P'F >i ) ]e i , 



e 



2 



\ J ' dXP"( V )F^e\ (2.13) 



0i and 62 represent the intrinsic and Sachs- Wolfe effects, respectively. The latter can be divided into the 
ordinary Sachs- Wolfe effect 6 sac and the Integral Sachs- Wolfe effect Qi sw , where 

&sac = \[{P"F) e -{P"F) l 

Q lsw = \f ' dXP"'{ v )F. (2.14) 
The second-order temperature perturbation is 

5T/T = / 1 (A e )[i/ 1 (A e )-r] - [A^' + t^] j\xA^ 

_ J^'dX^A^ + A^A^-A^' j\xA^\X)} + ^S^, (2.15) 
where (i], x l ) — (A, A — A) in the integrands and 

h(X e ) = ~ f dAP'i^eV, 

a {1) = - 2 i'i:,yr'. 

A ( V = -h [ PL{ + P 2 Mi + QN ij + Cf}e i e j . (2.16) 

Now let us assume an appropriate form of F as follows, to represent a situation including local inhomogeneities: 

F(x) = F P (x)+F L (x), (2.17) 

where the part of primordial density perturbations {F P ) and the part of local homogeneities (F L ) are expressed 
as 

F P = J dka(k)e 4kx , 

F L = R(r)Yr(e,<f>). (2.18) 

In the former equation, a(k) is a random variable, YJ m (#, (j>) is spherical harmonics, and (r,6,<p) is the polar 
coordinates, i.e. r 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . For the above F(x), we have the first-order density perturbation 

Sp/p = — 2 (-P' — l) (Ai^p + AF L ), 
1 pa z V a / 

1 to! 



(^ P ' - l ) [- J dka(k)fc 2 e lkx + £(r)Y; m (6», 0)] , (2.19) 



where 



with i? r = dR/dr. 



mn = ±l(r 2 R, r )- l «±VR ,2.20) 
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The part of AFp is for the ordinary primordial perturbations, so that by the averaging process, we have 

(a(k)a(k')) = (27r)- 2 7Mk),5(k + k'), (2.21) 

with 

Pjr(k) = PFofe" 3 (fc/*b) n_1 3?(fc), (2-22) 

where T s (fc) is the matter transfer function [15] and Vfo is the normalization constant. Then the average value 
of (<5i p/pf is 



((fp/p) 2 > 



par \ a 



(2tt)- 2 y dkP F {k)k A + {AF L 
~p^(i P ' ~ 1 )] 2 [( 2n y 2p Pof ^(fc/fc )"- 1 T s 2 (fc) + (AP L 



(2.23) 



The part of local inhomogencities represents a realization of cosmic variance (in our neighborhood) in primordial 
perturbations with small I. In order to consider directly the north-south asymmetry, we assume in the following 
a simplest case with I = 1 and m = 0, i.e. Fl — R(r) cos 9. As the observational background of this dependence, 
there is an asymmetric distribution of galaxies within 200 — 300 Mpc around us. According to the studies of 
galaxy number counts in the Sloan Digital Sky Survey (Yasuda et al.[16]), the galactic number density in the 
stripes toward the Southern Galactic Cap is larger than that in the stripes toward Northern Galactic Cap. Such 
an asymmetry of matter distribution may extend to the region on these scales. 
Then the first-order temperature perturbations are 



5T/T = 6p + e L cos0, 



(2.24) 



where 



6 P = 
Ql = 



- / dka(k) / d\P'(ri)(kn) 2 e 



2„ikx 



\ f ' d\P'( V )R r 



(2.25) 



and r in Eq.(2.11) was here neglected, because we pay attentions to the Sachs- Wolfe effect in the low-? cases 
such as I < 30. 

The former equation can be rewritten as 

Q P = f dk a (k){-i[(P'0 o + z^') OJ Pi(M)] + ^EH)'^ + 1 ) e ^) P '^)}' ( 2 - 26 ) 
J i 

where 

e P(I) = j * d\P'"(r ] )j l (kr) - {k(P') e (2l + l)" 1 ^ + l)j,+i(fcr e ) - lji-i(kr e )} + ji(kr e )(P") e }. (2.27) 

In these equations, we have r] = A and r = X — A. In the derivation of Eq.(2.27), we used the relations[17, 18] 

e *kx = e ikr„ = J2(-i) l (2l + l)j ; (fcr)P ; (/x) (2.28) 



and 



(21 + l)^Pi{p) = (l + l)Pi+i(p) + IPi-i(p), 



(2.29) 



where p = cos 9k and 9k is the angle between the vectors k % and e % . In the above equation for 0£, 9 is the angle 
between x 3 -axis and the vector e 4 , i.e. e 1 = sin 9 cos <f>, e 2 = sin sin 0, e 3 = cos9. 
The second-order temperature perturbation consists of three components: 



6T/T=e LP + e LL + &pp, 

2 



(2.30) 
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where ©lp is a component including the product of primordial perturbations and local inhomogeneitics, ©ll is 
a component including only the product of local inhomogeneities, and @pp is a component including only the 
product of primordial perturbations. In this paper we are concerned mainly with studying how CMB anisotropies 
are influenced by local inhomogeneities, and 6pp is small enough, compared with 6p in the first-order. From 
now, therefore, we neglect ©pp and treat only ©ll and ©lp- 

In Eq. (2.15) for S2T/T, the term with includes P^Mf^N^ and C{. The terms with and C{ are 

small, compared with the terms with P? and M|, because we have Q/P 2 < 10~ 2 always and the contribution 
of gravitational radiation is very small. So they are neglected in the following. 
For ©ll, we obtain from Eq. (2.15) 

©ll -©i°i + ©i 2 i cos 2 6, (2.31) 

where 

©l2 = ^ d\P'( V )r- 2 {R 2 + 2 -P[Vi(R r ) 2 + UR r R/r - 6(i?/r) 2 ]}, (2.32) 



q l1 = l(f e dXP'R,rr) 2 ~ \ f ' d\PP'{R, rr ) 2 + 1 I e d\P"R, rr I d\P{X)R. rr {\) 
+ - f K d\P'{--{R, r f RR, rr --(R/r) 2 

+ —P[7(R,rr) 2 - 7{R. r /r) 2 - Qr^R^R^ - R/r 2 ) - 6 r - 3 R(2R, r - i?/r)]}. (2.33) 
For ©lp, we obtain from Eq.(A23) in Appendix using the relation (2.29) and performing partial integrations 

Q LP /cos6=- f rfka(k)^(-i) ; (2/ + l)^ P P ; ( M ), 



(2.34) 



where 



n 



(0 

LP 



(P") e jl(kr e ) + [ ' dXP"'j l (kr)} x f ' d\P'(\)R ir (\) 
j * dXji(kr){p'P"R, rr + (P'") e PR, rr + ~(P'R, r ) 

+ P'R,rr + \^[ 1P ' R + PP'i^Rrr + 9i2, r /r - 9R/r 2 )} 



- -PP'k 2 (R^ rr - R, r /r + R/r 2 



d?_ 

dX 2 



P" J d\P(\)R,rr(X) } 



Ae 



dXP"(X)ji(kf) + I dXR, rr {kjl 1) (kr e )[(P') e P' + {P)eP"\ 



+ k 2 j\ 2 \kr e ){P') e P], 



(2.35) 



and we neglected the terms with <f>k in Eq.(A23) because they vanish in the k integration. Since the primordial 
random perturbations are rotationally symmetric and the local inhomogeneity is assumed to be proportional to 
cos#, it is reasonable that ©lp is represented as the product of cos# and rotationally symmetric perturbations 
(cf. Eq.(2.34). The latter can be expanded using the Legendre polynomials. 
In the above equations, we used the auxiliary functions defined by 



J i 1){kr) = zf—^ji-iW-il + l^ikr)}, 



f\kr) 



21 + 1 



(2l 2 + 2l-l)(2l + l) 1(1-1) , (j + m + 2) , 

Jiykr) - — -3i_ 2 {kr) —— Ji+2{kr) 



(2Z — 1)(2Z + 3) 



21 - 1 



21 + 3 



(2.36) 
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III. POWER SPECTRA OF CMB ANISOTROPIES 



The CMB anisotropies in the present analysis include the contributions from primordial perturbations and 
local homogeneities, and only the former perturbations are regarded as statistically random quantities. In order 
to derive the power spectra, therefore, we take the statistical average () only for the contribution from the 
primordial perturbations, and ((ST/T) 2 ) is expressed as 

((ST/T) 2 ) = ((ST/T) 2 ) + 2(ST/T5T/T), (3.1) 

1 1 2 

where we took first two terms and neglected higher order terms. For the first-order anisotropies, we have 

((5 T/Tf) = ((e P ) 2 ) + (Q L ) 2 cos 2 6, (3.2) 

where &l is defined in Eq.(2.25) and 

(To) 2 ((e P ) 2 )=^^±ic ; . (3.3) 

The power spectra C/ are 

Q = (T ) 2 J dkk 2 r F (k)\n^(k)\ 2 , (3.4) 

where To is the present CMB temperature and 

Hf(k) = -(P'% - k(P') e]l (kr e ) - (P") e]0 (kr e ) + d\P'"j (kr), 

n { p\k) = h[(P') -(P') e }j{ 1 \kr e )-(P") e j 1 (kr e ) + [ K d\P'".n(kr). (3.5) 

For / > 2, we have 

H P \k) = fc(P')ej i (1) (fcr e ) - (P") e ji(kr e ) + dXP^Mkr). (3.6) 



In the derivation of Ci, we used Eq.(2.21) for (a(k)a* (k')>, and the formulas for Pi(p) : J\ Pi(n)P v (n)dn = 
2/(2/ + 1),0 for 1 = 1', I ^ V, respectively. 
For the second term of Eq. (3.1), we have 

(ST/T ST/T) = (e P e LP ) +e L Q LL cos0, (3.7) 

1 2 

where 

e LL cos6» = e^cos0 + e^cos 3 e (3.8) 

and 

(e P e LP )/cos0 = i(2^)- 1 ^(2/ + i) f dkk 2 v F (k)nfnf L . (3.9) 

If we put (&p&lp) in the form 

(To) 2 (Opo L p) = Yl ^r ACi cos ( 3 - 10 ) 

i 



we have 



AC; = i(T ) 2 J dkk 2 V F (k)U P ) U% (3.11) 
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where Hp is expressed by performing the A differentiation in Eq.(2.35) as 

/■A e r\ e 



K» L = \ dX$ j l (kr) + {-(P") e j l (kr e ) + dXP"'ji{kr)]x dXP> (A)P r (A) 

A A A 

+ / ' d\P"R rr f dXP"(X)j l (kf) + f ' d\R rr {kf>(kr e )[(P') e P' + (P) e P"] 

+ k 2 3 \ 2) (kr e )(P') e P} 1 (3.12) 



where 



$ = - PP"P, 



2P"'p + Ip' - (P'") e P 



R 7 rr ~~ R,r — PR 



- \{WP" + PP'")(4Rrr + 9R.r/r-9R/r 2 ) + P"" ( d\PR. rr 

+ ^[{P'f + PP"](4P, rrr + 9P, rr /r - 18P r /r 2 + 18P/r 3 ) 

- -PP'(4R, rrrr + 9R, rrr /r - 27P jrr /r 2 + 54P ir /r 3 - 54P/r 4 ) 

+ yPP'/c 2 (P, rr -R, r /r + R/r 2 ). (3.13) 

It is found by numerical calculations that the dominant contribution to AC; comes from the last two terms 
(proportional to PP') in $, especially the terms with highest-order differentiations with respect to r. 

The total anisotropics in Eq.(3.1) consist of the rotationally symmetric component (((Op) 2 )), the asymmetric 
component ((0p0z,p)), and the dipole and quadrupolc components with (Ol) 2 and OlOll- The key point 
in this paper is that the north-south asymmetry affecting the CMB power spectra of I > 2 is caused by the 
component (0p6lp) which is proportional to cos#. 

IV. SIMPLE MODELS OF LARGE-SCALE LOCAL INHOMOGENEITIES 

In order to study the gravitational influence of local inhomogeneities on CMB, we consider a simple model 
with north-south asymmetry, in which Pl(x) is expressed as 

Pl(x) = R(r) cos 6. (4.1) 

This corresponds to the case I = 1 and m = in Eq.(2.18). For P, four types of functional forms are considered 
and compared : 

R= Po cxp[— a(x — l) 2 ], ip [l + cos2-7r(x — 1)], 

Rqx exp[— a(x — l) 2 ], and -P a; 2 [l + cos27r(a; — 1)] (4-2) 

in the interval x = [£1,2:2] with x = r/r c , where X\ = r\/r c — 0.5 and x 2 = r 2 /r c = 1.5. In all types we have 
P = for x > X2 or x < x\. Po is the normalization constant and a constant a is chosen as 20. Here a$r c is 
assumed to be 300/i _1 Mpc (Po = 100ft- km/s/Mpc). These types are called in the following the Gaussian 
type (G), the sine type (S), the modified Gaussian type (MG), and the modified sine type (MS), respectively. 
In the first two types, R is radially symmetric around the surface x — 1. In the second two types, P has small 
asymmetry around it. 

First, we derive the first-order CMB anisotropies for comparison. They are calculated using Eq.(3.4). In Tabic 
I, the behavior of 1(1 + 1)C/ is shown for n = 0.97. The root mean square (Ai) of 1(1 + 1)C;£ for I = 2 — 11 
is 0.172 for n = 0.97, where A 1 = {£;" 2 [^ + l)C{\ 2 /lQ} 1 ' 2 ^ £ = 27r/[Pp (T ) 2 ] and P is the present CMB 
temperature. The normalization constant Ppo can be determined as Ppo = [870(yuK) 2 /(2.7K) 2 ](27r)/^4i = 
2.1 x 10 -8 , so that it may be consistent with the observed CMB anisotropy. 

The first-order anisotropy 0^ due to local inhomogeneities is derived from Eq.(2.25) in the four types and 
the ratios of their values to Po in the four types are shown in Table II. 

Next we derive the second-order anisotropy for local inhomogeneities of the above four types, using Eq.(3.11). 
The behavior of AC/ is shown in Table I for I = 1 ~ 22. It is found that AC/ change the signs often in this 
interval of I , but their absolute values are comparable and increase slowly with the increase of I. The total 
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TABLE I: CMB anisotropy powers 1(1 + l)Ci and 1(1 + 1)AC; in the case n = 0.97. The latter is caused by the coupling 
of cosmological perturbations and local inhomogeneities of types G, S, MG and MS. Here £ = 2tt/[Pfo(To) 2 ], and Vfo 
and _Ro are the normalization factors. 

/(/ + l)Ci£ 1(T 3 x 21(1 + l)ACi£/Ro 
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TABLE II: CMB anisotropics caused by only local inhomogeneities of types G, S, MG and MS. Ro is the normalization 
factor. 



model types 


G 


S 


MG 


MS 


mean 


Ql/Ro 


0.75 


0.022 


0.92 


0.23 


0.48 


Qllo/(Ro) 2 


-9.1 x 10 3 


-9.1 x 10 2 


-9.1 x 10 2 


-8.7 x 10 2 


-9.0 x 10 2 


QlL2/(Rq) 2 


2.6 x 10 4 


1.6 x 10 4 


3.1 x 10 4 


2.4 x 10 4 


2.4 x 10 4 



power is C'i + 2AC; cos 9, and the north-south asymmetry is represented by the factor cos 9 in the second term, 
which is positive and negative for the northern and southern hemispheres, respectively. 

The ratios (A2/A1) are shown in Table III, where A 2 is defined as the mean square root (A2) of 21(1 + 1)AC;£ 
in I = 2 ~ 11, that is, A 2 = {J2]L 2 i 2l ( l + 1)AG] 2 /10} 1/2 ^. The dipole component (I = 1) is separately treated 
from the multiple components (I > 2). Here we remark that we have not specified the value of a^r c yet, and R 
is arbitrary. From Table III, we find that, as the mean value, we have 

\R \ = 7.5 x 10~ 5 (A 2 /A 1 ). (4.3) 
The second-order anisotropies Qllo/{Rq) 2 and @LL2/(-Ro) 2 are derived using Eqs.(2.32) and (2.33), and 



TABLE III: Ratios of second-order CMB anisotropies (A2) to first-order CMB anisotropies (Ai). The former is caused by 
the coupling of cosmological perturbations and local inhomogeneities of types G, S, MG and MS. Ro is the normalization 
factor. 



model types 


G 


S 


MG 


MS mean 


(A 2 /A 1 )/\R \ 


1.6 x 10 4 


0.81 x 10 4 


1.7 x 10 4 


1.2 x 10 4 1.3 x 10 4 
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their values in four types are shown in Table II. For the above values of Vfo and Ro, Ql, ©llo and <9ll2 are 
w 3.6 x 10- 5 (A 2 /A 1 ),-5.1 x IQ- 6 {A 2 /A 1 ) 2 and 1.2 x lO" 3 ^/^) 2 , respectively. 

In order to examine the significance of Rq given in Eq.(4.3), we consider here the first-order density pertur- 
bation due to local inhomogeneities. It is expressed as 

(6 p/p) L = \ i-P' -l)R cos 6, (4.4) 
l pa z \ a / 

where 

R = {r c )- 2 [(x 2 R, x ), x /x 2 - 2R/x 2 }. (4.5) 
The ratio (8M/M) of perturbed mass (5M) to background mass (M) in the interval x — [xi,x 2 ] is defined by 

f ' J (6plt>)Li?dxd»l{\l(x,f - (n) 3 ] / dp} 

The factor [(a'/a)P' — 1] can be regarded as having the value at present epoch, because the inhomogeneities are 
at the place of z = 0.05 ~ 0.15, and it is equal to —0.456 for the concordant background model with Q = 0.27 
and A = 0.73. Moreover, (pa 2 )o(r c ) 2 — 3[ar c /(c/H n )] 2 , where we have used the unit c = 8irG = 1 in this 
paper. Therefore we have 

J = -7.0(300/i _1 Mpc/a o r c ) 2 / [{x 2 R, x ), x - 2R]dx, (4.7) 

J X\ 

where x 2 — 1 = —{x\ — 1) = 0.5 and c/ H = 3000/i _1 Mpc. After performing the integration in Eq.(4.7), we 
obtain 

J = (2.4, 14.0, 13.6, 10.3) R a (SOO/i^Mpc/ao^) 2 , (4.8) 

for types G, S, MG and MS, respectively, and their mean value is J = 10.1 Ro (300/i~ 1 Mpc/aor c ) 2 . For i? 
given in Eq.(4.3), J leads to 

J = 7.6 x 10- 4 (A 2 /A 1 ) (SOO/i^Mpc/aor-c) 2 x R /\R \. (4.9) 

On the other hand, the observed value of (((5M/M) 2 ) 1 / 2 is about 1 on the scale of 8/i _1 Mpc and for the power 
spectrum P(k) cx k n , we have SM/M cx M _ (™+ 3 ^ 6 cx r -(™+ 3 )/ 2 = r -i.985 f or n _ q.97. Since we are considering 
a local inhomogeneity included in the sphere of radius 2aor c , its scale is regarded as 4ea r c (~ 1200e/i _1 Mpc) 
with e w 1, so that the value of SM/M for the inhomogeneity is 

(<5M/M) powor = (8fc- 1 /4eo r c ) 1 - 98B /& = 4.8 x 10- 5 (6e 1 - 985 )- 1 (300/i- 1 Mpc/aor c ) 1 - 986 , (4.10) 

where b is the biasing factor[19, 20]. If we assume \ J\ — (5M/M) powor , we obtain from Eqs.(4.9) and (4.10) 

Ail A x = 0.063(fee 1 - 985 )- 1 (300/i- 1 Mpc/a o r c ) 0015 . (4.11) 

For this value of A 2 /A\, we have 

|6 L | = 2.3x 10- 6 /(&e 1 - 985 ), e LLQ = -2.0xl0- 8 /(6e 1985 ) 2 , LL2 = 4.8x 10- 6 /(6e 1 ' 985 ) 2 for n = 0.97, (4.12) 

neglecting the factor of aor c . If foe 1 - 985 ^ 1 5 these are small, compared with \fC\/T (^ 5 x 10~ 5 ) and ^/C 2 /To(~ 
10 -5 ). If be 198,5 < 0.5, however, Qll2 is comparable with ^/C^/Tq or larger than it. 

Finally let us compare our above theoretical results with the observed anisotropy spectra by Eriksen et al. [4] . 
Here wc assume that the matter density perturbations in the Southern and Northern hemispheres are positive 
and negative, respectively, corresponding to the galactic number counts in SDSS[16]. Then we have Rq < in 
order that the directions 6 = and n may be in the Galactic North and South, respectively. In our results 
of calculations (from Table I), we find the trend that AC; for I = 2, 3 and 4 are positive for R < and AC; 
for / > 5 change their signs in the period of Al = 3. This trend seems to be seen in Fig. 2 of Eriksen et al.'s 
paper[4]. It is found, moreover, that corresponding to the north-south asymmetry (proportional to cos#) the 
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ratio ACi/Ci seems to be 0.3 ~ 0.5 by reading the points in Fig. 2 of their paper. From Eq.(4.11) and Eriksen 
et al.'s observational result, accordingly, we obtain the condition be 1985 = 0.1 <~ 0.2. This condition can be 
satisfied in the reasonable range of b. For instance we have e = 0.45 ~ 0.63 for b = 0.5. 



Under this condition, @ll2 is comparable with y/C^/To, so that the measured value of Ci may be disturbed 

by &LL2, because of their similar angular dependence. &LL2 and \J~C2~ depend on P 2 (cos#) and P2{cos9k), 
where 9 and 9k are the angles between a directional vector e and the x 3 -axis and between e and k, respectively 
(cf. Appendix), but both of them may be measured as the quadrupole components. So its measured value may 
have been given a value smaller than its theoretical expected value by the offset effect. 



In this paper we derived the asymmetry of CMB anisotropy powers (proportional to cos 9) with comparatively 
low I using the second-order perturbation theory by assuming the local matter distribution is dipole-like. This 
assumption seems to be rough, but it is consistent with the observed situation that the asymmetry of CMB 
anisotropics disappears when they are averaged in the whole sky and the asymmetry in the distribution of 
galaxies is also disappears similarly in the whole sky, though it has only a small north-south asymmetry. If we 
add a correction term oc cos 2 9 to Fl, more realistic simulations to observed AC; may be obtained, especially 
for I > 20. Moreover, to clarify the relation between the derived asymmetry and non-Gaussianity, it is necessary 
to investigate the multi-point correlations of anisotropies. Detailed analyses on these points are to be done in 
the next step. 

In our result, AC;/C; depend not so on the distance to the local inhomogeneity (i.e. the value of aor c ), but 
sensitively on be^ n+3 ^ 2 . It was found that the behavior of AC; in the case n = 0.90 is not so different from that 
in the case n = 0.97, so the conclusion in our result does not depend on n. 

The northern and southern hemispheres were assumed here as those in the Galactic coordinate frame. However 
our theory itself does not depend on any coordinate frames and can treat any local inhomogeneities which can 
be supposed. 

When the measurements of CMB anisotropies are more accurate, the relation between the local matter 
distribution and the CMB asymmetry will be more realistic through the second-order perturbation theory. 



The author thanks K.T. Inoue for helpful discussions on the asymmetry of CMB anisotropies. Numerical 
computation in this work was carried out at the Yukawa Institute Computer Facility. 
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APPENDIX A: DERIVATION OF 82 T/T 



From Eq.(2.15), we obtain 




(Al) 



neglecting the terms with r at the decoupling epoch. Each term in Eq.(Al) is expressed as follows: 




(A2) 



so that 



where the term (J'*" d\P' Fjje^i) is neglected by the reason described in the text. Similarly we obtain 




(A3) 



A^A^' = -PP'[(i^eV) 2 + 2F L , ij e i e> 



F PM e k e l ] 1 



(A4) 
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and 



where 



A«" J ' dXA^ = X -P" [F LM e k e l ' dXP(F Liij + F Ptl] )e 1 ^ + F PM e k e l ' ,1X1 ' F L ljl ' , , (A5) 

/ " dXA {2) ' = -\ [ ' dX[P'L{ + 2PP'M J i ]e i e j , (A6) 

Life, = —{Fl^) 2 - / /./' / ,,< V< + \f l ,iF l ,i - :!/'/.,<'//-..,<' 

- (F L F Pil] + F P F L<ij yei + ^F L<l F P , h 
Mle'e, = ^{[19F L , a F LJl - 12Ai^ L)ij ] e V - Z[F LM F LM ~ (AF L ) 2 } 

Zo 

+ [38F LiU F Ptjl - UAFlFp^j - 12AF P F L , y ]eV - 6[F LM F PM - AF L AF P }}. (A7) 
Using the radial coordinate r, we have 

F LJ e> = R r Yr, 

F PJ e j = J rfka(k)ifc^e ikx , (A8) 

corresponding to Eq.(2.18), where kx = kr/i, and 

F PJ F P,j = i J dka{\z)F LJ k 3 'e ikx , 
FpsFpji = - I dka(k)F L>jl k j k l e ik *. (A9) 



To express FljW and F^jik^k in terms of polar coordinates, we introduce two unit (three-dimensional) vectors 
e l g and together with e 4 , and make an orthonormal triad vectors with components: 

(e l ) — (sin cos 0, sin ^ sin 0, cos#), 
(e^) = (cos # cos <j(>, cos sin 0, — sin#), 

(e;) = (- sin 0, cos 0,0), (A10) 

respectively, where is the angle between the vector e % and the x 3 -axis. Using these vectors, we have dr/dx' 1 = 
Bi, d0/dx % — r~ 1 eei and d(j>/dx l — {rsva.6)~ 1 e ( j ) i for the polar coordinates (r,6,<j>). 

Moreover, the projection of k to the triad is expressed using another angles 9 k and <f> k as 

k l (eoi , e^i , e t ) = k (sin k cos 4> k , sin 9 k sin 4> k , cos 9 k ) , (All) 

where = ^e- 7 , etc, and the angle 0^ is the angle between the vectors e l and A; 1 . Then we have 

F L .tk l = F L j(e l k i e i + e l e k i e ei + e^k'e^) 

= k[F L ^ r cos9 k + F L fir~ x s\n9 k cos</> fc + F L ^(r sinO)' 1 su\0 k sin0 fc ], (A12) 

F L ,tjk l k j = k'k^FL^ej + {F Lfie + rF L , r )r 2 e w e 93 + {F L<H> + r sin 2 6F Lir 

+ sin9 cos9F Lj e)(r smO^e^e^ 

+ 2(F L . x e - r- 1 F L j)r~ 1 e i e 9j + 2(F L ^ 4> - r^ 1 F L ^)(rsmey 1 e i e 4 , j 

+ {F LM - cot 9F L ^){r 2 sin 0) _1 e w e#]. (A13) 

Substituting Fp and in Eq.(2.18) with i = 1 and m = into Eq.(A3), we obtain 

(/i) 2 = \(J ' dXP' R rr y cos 2 6 + 1 y " dXP' R rr cosO j ' dXP'i J dka(k)/c M e ikx , (A14) 

A^W' = ^PP' (i? rr ) 2 cos 2 6-2iR trr cosO J dka(k)(fc^)V kx , (A15) 
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J d\A^ = ip" fl.rry d\P(fj)R rr {f) cos 2 - R rr cos6 J dXP J dka{k){k^fe 

r-A 



2 ikx 



- y dka(k)(kfi) 2 e lWx J d\PR rr cos9 . (A16) 

For the terms in , we have 

FF lje l e j = PJ? rr cos 2 9 + J dka(k)e lkx [i? rr - [ky) 2 R] cos0, (A17) 

FjFi = (P, r ) 2 cos 2 + P 2 sin 2 6/r 2 + 2 y dka(k)e* kx ifc[P !r cos0^ - r _1 i? sin 6» sin 6» fc cos0 fc ], (A18) 

F, il F J ie i e j = (P, rr ) 2 cos 2 + (i? r ) 2 sin 2 0/r 2 + 2 j dka{k)e l ^[-R, rr cos9{k^) 2 

- R tr sin 9k 2 sin 9k cos(f>k}, (A19) 

AFF 13 eV = (P, rr + 2r~ 1 P r - 2r" 2 i?)i? jrr cos 2 

- j dka{k)e^[{R. rr + 2r- 1 R tr -2r- 2 R){kii) 2 + R^ rr k 2 ]cos9, (A20) 

FyFy - (AF) 2 = -4r" 1 ( J R ;T . - r - 1 R)(R yrr + r" 1 ^ - r~ 2 P) cos 2 + 2r- 2 (F, r - r^F) 2 

- 2y dka(k)fcV kx {P ! „. M 2 cos0 + r- 1 (P, r - r^i?) [cos sin 2 fc 

- 2sin0sin0 fe cos0 fc cos0 fe ] - (F :I . r + 2r~ 1 F, r - 2r~ 2 F) cos 0}. (A21) 

Using the above equations, we obtain 6^ and 6^p as follows: 

cA e ,o 1 /-A 



6ll = i(y C dAP'i? irr cos0) 2 -i y ' dXPP'(R. rr cos 9f 
+ \ f dXP"R, rr ( dXP(fj)R rr (f) cos 2 9 
+ 1 / * dAP' {-[##,„. + j|(F r ) 2 ] cos 2 + Jr- 2 i? 2 sin 2 



+ ^P[19(F,„.) 2 cos 2 + 19(F'/r) 2 sin 2 - UR. rr (R. rr + 2r~ 1 F, r - 2r~ 2 F) cos 2 
+ 12r- 1 (P, r -r- 1 P)(P, rr + r- 1 P lT . - r" 2 i?) cos 2 - 6(F, r - r^F) 2 ] | (A22) 



and 



n2 „ikx 



9 LP = ~\J " dXP'R.rr cos 9 x y ' dAP' y dka(k)(fc^i) 2 e 

+ iy ° dXPP'R. rr cos 9 y dka(k)(fc^) 2 e ikx 

- ^ 6 dXP" J dkaMikfife^^R.rr J d~XP(f])e lk5i + e lkx J dXP(fj)R, rr P(fj) 

+ ^ 6 dAP'{-^F, r cos0 i y dka(k)fc M e lkx - y dka(k)e 4kx [P !rr - (fc M ) 2 P]cos0 

- iy <2ka(k)e* kx r _1 Fsin0 iksin9 k cos0 fe | 

+ ^y rfAPP' J dka(k)e tkx {19[-R, rr cos 0(fc M ) 2 + F, r sin 0fc 2 /i sin fc cos<?i> fc ] 

+ 6fc 2 [(F rr + 2r _1 F ;T . - 2r~ 2 F)/i 2 + F >rr ] cos0 + 3fc 2 [F, rr cos 0^ 2 
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+ r 1 (i?, r — r 1 R)(cos 9 sin 2 9k — 2 sin 6 sin Ok cos Ok cos 4>k) — (R, rr + 1 R tT — 2r 2 R)cos9]} 



+ -(P')e / dka(k)(/c M ) 2 e ikx / d\PR, rr cosO. 



1 




(A23) 



By integrating Eq.(A23) partially with respect to A, we obtain Eq.(2.34). In the process of integrations, it is 
to be noticed that we use the boundary condition (R) e = (R) = and (d m R/dr m ) e = (d m R/dr m ) = for 
1 < m < 4. 
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